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ABSTRACT 


A theoretical  analysis  is  made  of  the  mechanical  effect  of  a detonation 
occurring  in  a closed  bomb  during  routine  propellant  testing  in  the  pressure 
range  for  which  convective  burning  is  expected  to  develop.  A hitherto 
unreported  phenomenon  was  found,  consisting  of  an  amplification  of  the 
peak  reflected  detonation  impulse  by  products  of  combustion  and  leading  to 
a maximum  impulse  as  a result  of  a detonation  occurring  during  the  first 
third  of  burning.  Graphical  results  of  the  analysis  are  presented  for  a 
standard  combustion  bomb  using  an  experimental  propellant.  Limitations 
imposed  by  instrumentation,  especially  pressure  gauges,  appear  to  be  more 
severe  than  these  imposed  by  mechanical  construction  over  the  range  of 
experimental  parameters  of  interest. 
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NOMENCLATURE 
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Y 
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Acoustical  velocity  in  gas 
Acoustical  velocity  in  bomb  material 


Isochoric  specific  heat  (propellant  and  fill  gas) 
Constants  in  detonation  equations 
Cj  - 3 x 454/(4nY)ln  30.48 

C2  * VAVr 
C3  * Vr  Vr 

C4  ' (tv.  - V/ro 

Cc  « y /p  R T 
5 o p u 0 

C,  B W.C  /R  T C 

O O VO  u o vp 

Specific  energy  of  reaction  of  test  propellant 
Weight  fraction  of  test  charge  burned 


Reflected  impulse  per  unit  area  in  Interaction  between  detonation 
wave  and  bomb  wall 

Constant  in  Eq.  (2)  for  resonant  period  of  bomb 

Mass 

Pressure 

Radius 

Universal  gas  constant  * 8.316  MPa  - co3/mol-*K 
Tensile  stresses  (definsd  in  text) 

Shear  stress 

Texipcraturc 

Time 

Velocity 

Voluue 


Yield  factor  to  give  TSf  equivalent 
Sealed  distance 
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Symbols, 

Y 

6 

y 

p 

°e.y 

°t,y 


Ratio  of  specific  heats 

Loading  density  of  propellant  into  bomb  5 

Molecular  weight 

Density 

Yield  stress,  shear 

Yield  stress,  tension  and  compression 

wall  thickness 


m /V. 
P o 


Subscripts 

a Ambient,  just  prior  to  detonation 

b Bomb 

c Combustion  conditions  (maximum  pressure) 

cr  Critical  conditions 

d Detonation  conditions 

e Effective 

f Expressed  in  feet 

i Inside  * 

lb  Expressed  in  pounds 

n Mean  value 

max  Maximum  value 

0 Outside  (radius);  with  no  combustion  (impulse);  otherwise  denotes 
fill  gas  conditions 

P Resonance  time 

p Propellant 

r Reference  conditions 

5 Scaled  value 

v IsocKoric 

1 Corresponding  to  maximum  I . 
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INTRODUCTION 


The  closed  combustion  bomb  is  regarded  as  a valuable  tool  for 
studying  deflagration  properties  of  solid  propellants  used  in  guns  and 
rockets.  It  is  desired  to  extend  the  pressure  range  of  such  studies  to 
provide  data  on  anomalous  combustion  resulting  either  from  changes  in 
the  basic  mechanism  governing  the  conductive  regression  rate,  or  from 
development  of  burning  in  physical  grain  defects  - the  so-called  convec- 
tive burning  mode.  Tests  would  be  carried  out  over  a range. of  initial 
pressures  and  loading  densities  in  order  to  achieve  control  over  pres- 
sure range  and  pressurization  rate.  There  is  no  assurance  that  anomalous 
burning,  considered  as  the  first  stage  of  transition  from  deflagration 
to  detonation  (DDT),  can  be  limited  to  deflagration  under  test  condi- 
tions. Therefore,  it  is  necessary  Co  consider  the  design  criteria  for 
containing  a detonation  if  it  should  occur  during  a routine  deflagration 
study.  This  report  presents  principles  and  an  illustrative  example 
related  to  determination  of  requirements  for  containing  peak  deflagra- 
tion pressures  and  impulsive  detonat-on  loadings  for  elosed  bombs  over  a 
range  of  design  and  operating  parameters.  The  study  is  limited  to 
spherical  bombs,  for  which  the  theory  is  better  established  than  for 
other  shapes.  Results  are  based  on  typical  data  and  arc  intended  to 
indicate  trends  rather  than  establish  a final  design,  which  would  require 
additional  calculations  relevant  to  actual  construction  materials  and 

propellants. 

♦ 

The  purpose  of  this  study  is  to  present  conservative  design  princi- 
ples in  order  to  avoid  destroying  or  damaging  material  facilities. 

Unde r no  circumstances  eficutd  caution  he  lessened  uhen  /tutsan  safely  is 
involved;  remote  opeiu lions  are  aluays  advised  v/sen  high  explosives  cure 
concerned. 


DESIGN  PRINCIPLES 


Typical  schematic  pressure-time  histories  for  detonation  and  defla- 
gration in  a closed  bomb  are  shown  in  Figure  1.  The  pressure  at  time  A 
is  tins  initial  pro-pressure  provided  for  the  test.  Normal  deflagration 
proceeds  along  the  line  ABDE,  the  pressure  rising  in  a nearly  spatially 
uniform  manner  from  pA  to  pE>  The  portion  of  the  curve  BCD  indicates 
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FIGURE  1.  Schematic  Pressure-Time  History 
for  Closed  Bomb  Detonation  and  Deflagration. 


the  effect  of  a detonation  occurring  at  tine  tg  after  an  initial  period 
of  normal  burning  which  may  even  approach  zero.  Because  of  the  propaga- 
tive nature  of  the  detonation,  the  pressure-time  history  shown  as  BCD 
would  vary,  depending  on  the  radial  distance  from  the  detonating  source. 
It  is  convenient  to  consider  three  regions  surrounding  the  detonation; 
ur  decorations  at  standard  atmospheric  conditions  these  are:  (1)  the 
region  of  direct  explosion  effects,  extending  from  the  center  of  the 
explosion  to  a radial  distance  of  1.6  charge  radii,  and  characterised  by 
peak  pressures  varying  from  several  hundred  kilobars  down  to  several 
hundred  bars;  (2)  the  region  of  composite  blast  (Involving  products  and 
ambient  gas),  extending  outward  to  a distance  of  approximately  16  charge 
radii  and  peak  pressure  of  173  psia;  and  (3)  the  region  of  air  shock 
beyond  16  charge  radii  in  which  the  shock  pressure  eventually  decays 
with  increasing' distance  to  a simple  sound  wave.  It  should  be  noted 
that  in  an  unconstrained  environment,  the  line  A£  would  be  horizontal. 

A tractable  analysis  is  attainable  only  for  simple  geometries.  The 
simplest  is  the  sphere  with  a centrally  situated  spherical  charge  so 
that  the  expanding  shock  reaches  all  portions  of  the  wall  simultaneously 
without  the  shock  reinforcement  afforded  by  angular  configurations. 
Departures  from  this  ideal  configuration  have  not  been  treated  compre- 
hensively by  any  quantitative  analysis  and  would  result  in  more  severe 
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requirements  for  containment;  it  is  estimated  that  at  an  additional 
safety  factor  of  1.5  would  be  needed  for  a cylindrical  bomb  with  an  L/D 
ratio  of  unity. 

In  order  to  avoid  the  extremely  high  explosion  pressures  of  region 
one,  loading  densities  are  limited  to  values  below  0.A  g/cm^  and  associ- 
ated static  pressure  rises  of  approximately  90  kpsi.  The  final  pressure 
Pg'is  readily  calculated  from  thermochemical  and  gas  law  considerations; 
standard  uhfired  pressure  vessel  design  principles  suffice  for  estimat- 
ing static  pressure  mechanical  requirements.  Dynamic  requirements 
depend  on  determination  and  analysis  of  the  curve  ABODE  appropriate  to 
the  bomb  radius.  The  two  following  sections  of  this  report  will  con- 
sider the  dynamic  and  long  range  static  aspects  of  detonation  in  a 
closed  bomb.  Dynamic  or  detonation  effects  will  be  discussed  first 
because  the  static  pressure  conditions,  existing  after  decay  of  the 
detonation  transient,  can  be  obtained  as  a limiting  case  of  one  of  the 
equations  required  to  determine  detonation  conditions.  Except  for 
different  heat  loss  effects  and  possible  non-equilibrium  chemistry 
accompanying  detonation,  the  final  static  pressure  is  independent  of  the 
occurrence  of  detonation. 


DETONATION  EFFECTS 

While  internal  pressure  alone  dictates  the  design  of  a closed  bomb 
to  contain  static  loads,  the  situation  is  entirely  different  for  the 
dynamic  loading  presented  to  the  bomb  wall  by  a detonation  shock.  It 
has  been  found  that  the  effect  of  a shock  wave  depends  primarily  upon 
the  impulse  per  unit  area  (integral  of  the  pressure-time  curve)  and  its 
time  of  delivery  in  the  interaction  between  shock  and  wall.  When  a 
critical  Impulse  is  delivered  to  the  bomb  wall  within  a critical  time, 
it  may  be  expected  that  the  bomb  will  suffer  unacceptable  damage.  Four 
quantities  must  therefore  be  determined  in  order  to  estimate  the  response 
of  the  bomb  to  an  internal  detonation:  (1)  actual  impulse,  (2)  actual 
action  time,  (3)  critical  Impulse,  and  (A)  critical  action  time.  Items 
(1)  and  (2)  represent  detonation  parameters  which  depend  on  the  size  and 
energy  of  the  detonating  charge,  the  distance  of  the  bomb  wall  from  the 
charge,  and  the  nature  of  the  intervening  medium  just  prior  to  detona- 
tion. The  critical  parameters,  items  (3)  and  (A),  depend  upon  material 
properties,  especially  under  conditions  of  dynamic  loading,  and  upon 
bomb  geometry. 

Detonation  Parameters 


Items  (1)  and  (2),  mentioned  above,  are  the  detonation  parameters 
which  arc  needed  tb  determine  the  ability  of  a closed  bomb  to  withstand 
an  internal  detonation.  A third  quantity,  the  peak  reflected  pressure, 
will  also  be  calculated  as  an  indication  of  the  requirement  of  instru- 
mentation to  withstand  short  duration  pulses  of  extremely  (iigh  magni- 
tudes compared  to  the  static  pressure  capabilities. 
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Thd  bases  of  calculation  of  detonation  parameters  are  model  experi- 
ments and  scaling  laws  derived  from  dimensional  analysis.  Figure  2, 
adapted  from  Cohen,  shows  typical  relationships  between  scaled  para- 
meters pertinent  to  the  problem  at  hand.1 2 * 4  A discussion  of  the  rationale 
of  scaling  laws  is  contained  in  reports  by  Kinney  and  Baker,  while  the 
use  of  alternate  sets  of  parameters  (leading  to  similar  results)  is 
covered  by  Rinehart  and  Pearson.2-4  It  is  important  to  note  that, 
usually,  scaled  parameters  are  not  dimensionless  groups,  but,  as  in 
Figure  2,  have  dimensions  determined  by  the  units  in  which  the  model 
experiment  was  conducted.  While  it  would  be  possible  to  convert  the 
number  scales  on  the  coordinate  axes  to  agree  with  the  SI,  c.g.s.,  or 
some  other  standard  system  rather  than  the  English,  it  is  more  conven- 
ient to  convert  the  quantities  with  which  the  figure  is  entered  from 
c.g.s.  to  English. 

The  scaled  parameters  of  Figure  2 contain  four  scale  factors  (see 
the  Appendix  for  details):  acoustic  velocity  ratio  aa/ar>  density  ratio 
oa/cr,  pressure  ratio  pa/pr,  and  charge  weight  nu.  The  last  is  also  a 
ratio,  but  since  the  refeience  weight  in  Figure  2 is  1 lb  of  TNT,  the 
ratio  is  numerically  equal  to  the  actual  weight  in  pounds. 

It  1«  convenient  for  computational  purposes  to  express  the  scaled 
parameters  in  terms  of  design  parameters  of  the  bomb,  namely,  radius, 
loading  density,  and  fraction  of  propellant  mass  burned,  as  well  as 
other  constants  characterizing  the  initial  fill  gas  and  propellant 
combustion  products.  An  inert  fill  gas,  such  as  nitrogen,  will  be 
assumed,  thereby  eliminating  additional  chemical  heating  from  Interac- 
tions of  fill  gas  and  propellant  products.  Details  of  this  transforma- 
tion are  shown  in  the  Appendix.  The  ambient  conditions,  denoted  by  sub- 
script a,  refer  to  the  conditions  of  the  bomb  contents  at  the  instant  of 
detonation.  They  are  a function  not  only  of  the  fill  gas  parameters, 
but  of  the  fraction  of  the  propellant  charge  burned  at  the  time  of 


lEdward  Cohen  (consulting  cd.).  "Prevention  and  Protection  Against 
Accidental  Explosion  of  Munitions,  Fuels,  and  Other  Hazardous  Mixtures," 
Annals  of  the  Hen  fork  Academy  of  Sciences,  Vol.  1S2,  Art,  1,  p.  320. 

2Naval  Weapons  Center,  Engineering  Elements  of  Explosions,  by 
G.  F.  Kinney,  Naval  Postgraduate  School.  China  Lake,  Calif,,  NWC, 
November  1968.  (NWC  TP  4656,  publication  UNCLASSIFIED*.) 

University  of  Texas  Press.  Exp  ’oeions  in  Air , by  Wilfred  E. 

Baker.  <.ustln,  Texas,  1973. 

4 American  Society  for  Metals.  Sehaoior  of  Metals  Under  Itxpulsive 
Loads,  by  John  S.  Rinehart  and  John  Pearson.  Cleveland,  Ohio,  1954. 
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detonation.  This  latter  effect  distinguishes  the  current  study  from 
available  calculations  neglecting  detonation  augmentation  resulting  from 
partial  combustion. 

A computer  program,  incorporating  Figure  2 and  the  expressions 
developed  in  the  Appendix,  was  used  to  calculate  Ij,  p^,  and  t<j  as  func- 
tions of  f(0-l),  6(0.005-0. 3g/cm3) , and  po(0.1-30  MPa).  Propellant  gas 
properties  were  obtained  by  running  the  NWC  propellant  evaluation  pro- 
gram (PEP)  for  ALTU-16,  an  experimental  propellant  considered  for  the 
actual  tests;  the  fill  gas  was  assumed  to  be  nitrogen;  and  the  reference 
conditions  applicable  to  Figure  2 are  air  at  one  atmosphere  and  70°F 
(294  K).  A bomb  ID  of  5 inches  was  used,  corresponding  to  a commercial- 
ly available  bomb.  A summary  of  the  constants  used  is  provided  in  Table 
1.  As  mentioned  before,  impulse  is  one  of  two  Important  detonation 
parameters;  a noteworthy  result  of  the  calculations  is  the  existence  of 
a maximum  impulse  in  the  Ij  vs  f relationship  as  a result  of  competing 
effects  of  charge  weight  and  ambient  gas  conditions.  This  effect  prt  * 
vails  over  the  range  of  parameters  investigated  and  is  shown,  for  illus- 
trative purposes,  in  Figure  3 for  specified  values  of  loading  density 
(6)  and  fill  gas  pressure  (p0).  The  maximum  and  corresponding  value 
of  f vary  with  6 and  p0;  maximum  1^  is  show  as  a function  of  p0  and  6 
in  Figure  4 over  the  range  of  values  investigated.  Also  indicated  (for 
6 « 0.005  and  0.3)  are  typical  impulse  values  calculated  by  assuming  no 
combustion,  i.e.,  the  entire  charge  detonates  at  the  Initial  fill  gas 
pressure.  It  is  seen  that  the  trends  as  a function  of  fill  gas  pressure 
are  opposite  for  the  two  methods  of  calculation.  If  impulse  from  a 
detonation  under  initial  conditions  (no  combustion,  f«0)  is  compared 
with  the  maxiBHim  impulse  (f  ranging  from  0-0.3  depending  on  parameters), 
an  impulse  amplification  factor  (Idl^do>  observed  which  ranges  up  to 
nine.  Figure  5 is  a plot  of  this  factor  as  a function  of  6 and  p0  over 
the  range  of  parameters  studied.  Since  the  instant  of  DDT  is  not  yet 
predictable,  one  should  base  the  detonation  capabilities  of  the  boob  on 
maximum  impulse. 

The  second  ioportant  detonation  parameter  is  the  time  duration  of 
the  reflected  impulse,  shown  in  Figure  6 as  a function  of  loading  den- 
sity and  fill  gas  pressure.  The  relevance  of  tj  will  become  clear  in  a 
later  section. 

Although  peak  detonation  pressure  is  not  included  among  the  design 
requirements  for  the  bomb,  it  is  important  in  assessing  requirements  and 
capabilities  for  instrumentation,  especially  pressure  transducers.  Like 
impulse,  detonation  pressure  has  a maximum  when  considered  as  a function 
of  f;  however,  maximum  impulse  usually  results  froo  a detonation  which 
occurs  earlier  (f*0-0. 3)  in  burning  than  does  maximum  detonation  pressure 
(f*0. 6-0.7).  Figure  7 presents  maximum  calculated  detonation  pressures 
vs.  initial  fill  pressure  and  propellant  loading  density  over  the  range 
of  parameters  of  this  study. 
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TABLE  1.  Physical  Properties  and  Conditions 
Used  in  Calculations. 


Reference  conditions  (air) 


Pressure,  p 
Temperature,  Tr 
Molecular  weight, 

Initial  conditions 


Pressure,  p 

o 

Temperature,  T0 

Molecular  weight  (N^) , 

Specific  heat,  C ° 

r vo 

Propellant  gas  (ALTU-16) 

Isochoric  flame  temperature,  T 
Molecular  weight  ^ 

Specific  heat,  Cv 

Solid  propellant 


Density,  pp 
Yield  factor 


V 


0.10136  MPa 
294  K (70°F) 
29 


Variable 
240  K (70°F) 

28 

0.1786  cal/g  K 


4295  K 
19.73 

0.503  cal/g  K 


1.75  g/cm^ 


1.3 


The  ambient  pressure  at  all  burnt  (fnl)  is  the  pressure  which  would 
be  generated  in  the  absence  of  detonation  and  would  also  be  the  pressure 
after  detonation  transients  decayed,  assuming  similar  heat  losses  and 
products  of  reaction  in  the  two  cases.  We  consider  this  pressure  now 
because  it  is  the  natural  limit  of  no  detonation.  Figure  8 shows  the 
static,  long  term,  combustion  pressure  as  a function  of  6 and  pD.  The 
effect  of  p0  is  mainly  due  to  dissimilarities  in  properties  between  the 
fill  gas  and  pfopellant  combustion  products.  The  combustion  pressure  is 
used  to  determine  static  design  as  will  be  presented  in  a subsequent 
section. 

This  concludes  the  discussion  on  calculation  of  detonation  param- 
eters. The  nert  section  considers  the  determination  of  critical  values 
of  impulse  and  time,  necessary  to  the  estimation  of  the  capabilities  of 
the  bomb  under  detonation  loading. 
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FIGURE  3.  Relation  of  Reflected  Impulse  in  a Closed  Bomb  to 
Fraction *of  Fropollant  Burned,  Loading  Density,  and  Fill  Pres- 
sure. (ALTU-16  propellant,  5-inch  ID  spherical  bomb.) 
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FIGURE  6.  Impulse  Duration  Corresponding  to  Maximum  Impulse  for  Detona- 
tion in  Closed  Boob.  (AlTU-16  propellant,  5-inch  ID  spherical  bomb.) 


FIGURE  8.  Maximum  Combustion  Pressure  in  a 
Closed  Bomb  (ALTU-16  Propellant). 


Critical  Parameters 


In  order  for  a detonation  to  cause  failure,  the  reflected  impulse 
and  its  time  duration  must  bear  a certain  relationship  to  critical 
values  of  Impulse  and  time  which  depend  upon  material  and  structural 
details.  The  criterion  is:  impulse  must  exceed  the  critical  value  and 
must  be  delivered  in  a time  not  exceeding  the  critical  time.  The  prin- 
ciples of  calculation  of  the  critical  values  are  presented  in  the  fol- 
lowing sections. 

Critical  Impulse.  The  critical  impulse  is  related  to  the  critical 
particle  velocity  of  the  resisting  structure,  which  is  actually  that 
velocity  to  which  one  portion  of  a structure  must  be  accelerated  rel- 
ative to  another  in  order  that  failure  may  occur.  This .velocity4  (Ucr) 
for  steel  is  around  200  ft/s  and  Is  related  to  critical  impulse5  by  the 
expression 


*Naval  Weapons  Center.  Response  of  Structures  to  Blast:  A New 
Criterion , by  Robert  G.  S.  Sewell  and  G.  F.  Kinney.  China  Lake,  Calif., 
NWC,  June  1968.  (NWC  TP  4422,  publication  UNCLASSIFIED*) 
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3 

where  pj,  is  the  density  of  the  metal  in  slugs/ft  and  x is  the  thickness 
of  the  resisting  medium  in  feet.  For  steel,  a good  value  for  Icr  is 
1750  psi-ms/in.,  which  means  that  an  impulse  of  1750  psi-ms,  applied  to 
a 1-inch-thick  steel  plate,  will  cause  failure,  jLf  the  time  duration  of 
the  pulse  does  not  exceed  a critical  time. 


Critical  Time.  The  concept  of  critical  impulse  within  a critical 
time  is  discussed  by  Kinney  and  Sewell. 2 Both  the  results  of  detailed 
analysis  and  physical  reasoning  suggest  a critical-  time  equal  to  one- 
quarter  of  the  natural  resonance  period  of  the  structure.  The  physical 
argument  is  that  maximum  displacement  occurs  if  the  impulse  is  delivered 
during  the  interval  of  increasing  displacement,  i.e.,  the  first  quarter 
period.  For  a thick  spherical  shell,  loaded  symmetrically,  the  funda- 
mental frequency  is  given  by  Seide.6  The  corresponding  quarter  period 
is  given  by 


(2) 


where  \ is  the  mean  radius,  C is  the  velocity  of  sound  in  the  struc- 
tural material,  and  k is  a constant  depending  on  Poisson's  ratio  and 
thickness.  Considering  the  accuracy  of  the  other  available  data,  it  is 
sufficient  to  use  k a 0.75.  The  calculated  quarter  periods  will  then  be 
in  error  by  no  more  than  5X  for  steel  with  a Poisson  ratio  of  0.33 
regardless  of  wall  thickness. 

Knowing  the  critical  impulse  and  time,  we  can  now  determine  the 
effective  impulse.  If  impulse  action  time,  t({,  is  less  than  or  equal  to 
a quarter  period,  all  the  impulse  is  delivered  within  the  required 
critical  time  and  the  effective  impulse  is  equal  to  the  actual  Impulse. 
If  the  action  time  is  greater  than  a quarter  period,  only  a portion  of 
the  actual  impulse  is  delivered  within  the  required  critical  time.  The 
calculation  of  effective  impulse  is  then  determined  by  assuming  the 
pressure  time  curve  to  decay  linearly  (adequate  accuracy  compared  to 
other  assumptions)  so  that 


**Pnul  Selde.  "Radial  Vibrations  of  Spherical  Sheila,"  Journal  of 
Applied  mVa,  37,  528-530  (1970). 
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Having  now  completed  the  discussion  of  detonation  principles,  we 
turn  next  to  the  relatively  simpler  matter  of  requirements  imposed  by 
static  pressure  in  a spherical  pressure  vessel. 


STATIC  PRESSURE  CAPACITY  OF  THICK  WALLED  SPHERES 


The  pressure  capability  of  a thick  walled  sphere  is  presented  in 
numerous  references  based  on  various  criteria  which  are  still  subject  to 
controversy.  The  thick  walled  vessel  is  complicated  by  the  non-uniform 
distribution  of  stresses  across  the  wall  and  associated  combined  stresses. 
We  have  chosen  the  following  formulas  as  a basis -for  design:7 


8 

max^ 


s 

toax2 


a P 


3 3 

Rq  + 2R1J 

2<R03  - Rj3) 


s 

oax3 


s ° p 
s max  v 


(6) 

<7) 


The  quantities  sj,  82 « and  S3  arc  the  meridional  wall  stress,  hoop  wall 
stress,  and  radial  wall  stress,  while  ss  is  the  shear  stress.  All  four 
stresses  are  maximum  at  the  inner  radius,  Rj,  Important  parameters  are 
the  working  strength  in  tension  (assumed  equal  to  compression)  and  shear 
(usually  taken  to  be  0.6  of  tensile).  Based  on  the  foregoing,  shear 
stress  is  limiting  for  Rj/Ro  less  than  0.91,  while  tensile  stress  is 
limiting  for  Ri/Rq  greater  than  0.91.  Since  Ri/R0  is  less  than  0.91  for 
most  thick  wall  vessels,  the  static  pressure  design  will  be  based  on  Eq. 
(7)  with  Sjj  max  equated  to  the  working  shear  stress  for  the  material  in 
question. 


This  completes  the  section  on  design  principles  for  containing 
impulsive  and  static  loads  resulting  from  a detonation  occurring  Inside 
a spherical  bomb.  We  now  proceed  to  the  application  section  where  we 
will  examine  the  capabilities  of  a given  bomb  which  would  be  suitable  as 
a vehicle  for  studying  development  of  convective  burning  at  moderate 
pressures. 


* 


. 7McCrav-Hill  Book  Company.  Formulas  For  Stress  and  Strain,  by 
Raymond  J.  Roark.  New  York,  N.Y.,  1965. 
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APPLICATION 


The  purpose  of  this  section  is  to  illustrate  the  principles  just 
developed  by  a practical  example  in  which  the  operating  limits  are  to  be 
estimated  for  an  available  high  pressure  bomb,  whose  specifications  are: 

Inside  diameter  - 5 inches 

Outside  diameter  - 10  inches 

Material  - 316  stainless  steel 

Working  Pressure  - 16,000  psi 

It  is  desired  to  know  the  allowable  limits  of  propellant  loading  density 
to  keep  within  design  specifications  for  static  pressure  as  well  as 
impulsive  loading.  Also,  it  would  be  interesting  to  know,  as  a check- 
point, if  the  specified  working  pressure  is  over-conservative.  The 
tensile  yield  stress  for  316  stainless  steel  is  found  to  be  50,000  psi.9 
Using  os>y/ot  y = 0.6  and  a safety  factor  of  2,  we  calculate  a working 
shear  strengtn  of  15,000  psi.  From  Eq.  (7),  the  working  pressure  is 
calculated  to  be  17,500  psi  so  that  the  16,000  psia  specified  working 
pressure  is  reasonable. 

The  critical  impulse,  at  1750  psi-ms/in.  of  wall  thickness  is 

I - 2.5  x 1750  - 4375  psi-os  (8) 

The  critical  time,  given  by  Eq.  (2),  using  15,000  ft/s  as  the  acoustic 
velocity  in  steel,  is 

t « kR  /C  - 0.75(3. 75/12)/15000  « 15.6  us  (9) 

cr  in 

Operational  limits  are  determined  from  the  working  pressure,  criti- 
cal impulse,  and  critical  time  just  computed.  The  principles  involved 
are: 

1.  Working  pressure  is  the  limiting  consideration  for  combina- 
tions of  6 and  p0  which  produce  a final  static  pressure  greater 
than  16,000  psi  (110  MPa)  as  determined  from  Figure  8. 

4 

2.  Impulse  is  limiting  for  combinations  of  6 and  pD  which  produce 
a maximum  impulse  as  determined  from  Figure  4 which  leads  to 

an  effective  impulse  of  4375  psi-ms.  In  this  example.  Figures  4 
and  6 show  that  whenever  impulse  exceeds  4375  psi-ms,  the 


eAir  Force  Materials  Laboratory.  Aerosmce  Structural  Metals  Hand- 
book,. Vol,  1,  Ferrous  Alloys.  Wright-Pattersoa  AFB,  Ohio,  1964. 
(Publication  UNCLASSIFIED.) 
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action  time  is  less  than  the  critical  time  of  15.6  us,  so  that 
impulses  are  always  1002  effective  when  they  are  limiting. 

Limiting  combinations  of  6 and  p0  are  shown  for  static  and  dynamic 
criteria  in  Figure  9.  The  region  below  the  line  ABC  is  available  for 
conducting  experiment  ■}.  It  is  seen  that  static  considerations  are 
limiting  at  high  fill  pressures,  while  dynamic  considerations  are  limit- 
ing at  low  fill  pressures.  Lines  of  constant  maximum  static  and  detona- 
tion pressures  are  also  indicated  in  the  region  which  is  open  to  experi- 
ment. 


FIGURE  9,  Performance  Capabilities  of  5-Xncb  ID 
Spherical  Bomb  Using  ALTU-16  Propellant. 
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Two  additional  curves  are  shown,  labeled  "attainable  static  limit" 
and  "R0  =5.5  inches."  The  first  is  the  limiting  static  performance 
capability  of  a 5-inch  ID  sphere,  regardless  of  the  wall  thicknesss,  and 
corresponds  to  the  limiting  fora  taken  by  Eq.  (7)  when  Rq  becomes  very 
large.  The  second  curve  shows  the  improvement  in  dynamic  performance 
afforded  by  increasing  the  wall  thickness  from  2.5  to  3.0  inches. 

Figure  9 is  a summary  of  results  as  well  as  a typical  design  chart 
giving  relations  between  experimentally  chosen  parameters  and  per- 
formance. The  limit  line  ABC  not  only  delineates  the  region  of  avail- 
able combinations  of  pQ  and  6 for  safe  operation,  but  separates  propel- 
lant dependent  quantities  from  bomb  dependent  quantities.  The  bomb 
design  does  not  influence  the  lines  of  constant  maximum  reflected  detona- 
tion pressure  (pj)  and  constant  maximum  combustion  pressure  (pc)  which 
depend  solely  on  fill  gas  and  propellant  gas  properties.  Likewise  the 
result  of  changing  the  bomb  design  would  only  change  the  limit  of  the 
region  of  safe  operation,  i.e.,  the  line  ABC.  The  dynamic  limit  line  AB 
could  theoretically  be  increased  indefinitely  by  thickening  the  wall; 
improvement  in  the  static  limit  is  constrained  by  Eq.  (7)  unless  con- 
struction materials  of  higher  shear  strengths  arc  used.  However,  the 
limiting  feature  in  the  currently  available  5-inch  ID  design  is  more 
likely  the  high  detonation  pressure  peaks,  even  though  they  have  a 
duration  of  only  n few  microseconds.  It  will  probably  be  necessary  to 
interpose  a low  pass  filter,  cutting  off  below  100  kHz  in  order  to 
protect  any  sensitive  pressure  transducer.  If  this  can  be  accomplished, 
then  a loading  density  of  0.02  g/cm3  (20-g  charge)  will  provide  a 
maximum  test  pressure  ranging  between  35  and  100  MPa  (5075  and  14,500 
psia)  depending  on  initial  pressurization.  Although  the  propellant  data 
leading  to  Figure  9 are  based  on  a particular  propellant,  it  would  be 
simple  to  construct  analogous  graphs  for  other  propellants. 
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Appendix 

TRANSFORMATION  OF  SCALED  DETONATION  PARAMETERS 
TO  BOMB  PARAMETERS 


The  correlations  of  detonation  parameters  upon  which  Figure  2 is 
based 1 were  originally  determined  for  ambient  conditions  corresponding 
to  air  at  sea  level  pressure  and  70°F  (294  K) . In  order  to  apply  this 
information  to  the  current  study,  two  transformations  were  required. 
First,  a scaling  to  convert  to  different  ambient  conditions  (determined 
by  partial  combustion  of  the  propellant  charge)  was  applied  in  accordance 
with  the  discussions  in  Ref.  2.  Second,  the  resulting  scaling  laws  were 
expressed  in  terms  of  bomb  parameters  and  properties  of  the  propellant 
and  fill  gases.  The  result  of  the  first  transformation  is  shown,  without 
discussion  (see  Ref.  2),  in  the  following  equations. 
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(A-4) 


Subscripts  on  the  symbols  for  weight  and  radius  emphasize  the  dimensional 
nature  of  the  scaled  parameters.  In  addition,  the  subscript  c,  on 
charge  weight,  implies  that,  for  the  propellant  considered,  the  equiva- 
lent weight  of  the  reference  explosive  TNT  should  be  used.  This  point 
will  be  covered  later. 


' Calculation  of  the  actual  detonation  parameters  (subscript  d)  from 
the  scaled  values  (subscript,  s)  require  the  calculation  of  ratios  of 
ambient  to  reference  conditions  for  pressure,  density,  and  acoustic 
velocity  in  teres  of  experimental  parameters.  In  addition,  changes  in 
the  system  of  units  from  'English  to  metric  are  applied  to  scaled  distance 
and  pressure. 

The  first  factor  , in  the  scaled  distance,  is  converted  to 

metric  by  using  1 AD,€ 
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- Rf  = R1/30.48  (A-5) 

m,.  - 4/3  6(1  - f)Y/454  (A-6) 

lb,e  l 

where  Y is  the  yield  factor  for  converting  actual  propellant  weight  to 
TNT  equivalent.  We  have  as  a result 

V°lb!e  * C1/61/3(l  - f)1/3  (A-7) 


The  density  ratio,  pa/pr,  which  occurs  repeatedly  as  a scale  factor, 
is  expressed  in  terms  of  experimental  parameters,  thus: 


m + f6Vb 

Vb  - « - f)SVb/pp 

Vo  Vb 
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and 
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Combining  Eq.  (A- 8)  through  (A-10)  provides  the  desired  expression 


C2  *o<1  * 4 V * V* 

Pr  1 - (1  - f )6/pp 


(A-ll) 


The  perfect  gas  law  is  used  throughout  the  analysis  as  an  adequate 
approximation  of  the  equation  of  state  for  the  pressures  involved  in 
view  of  the  other  assumptions  which  are  made. 

The  pressure  ratio  is  most  easily  determined  indirectly  in  terms  of 
temperature  and  molecular  weight  ratios.  For  temperature  ratio,  Ta/Tr, 
the  ambient  temperature  is  deduced  from  the  energy  equation,  using  the 
conservative  assumption  of  no  heat  loss. 

fm  Ac  ♦ <t  - T )(fm  C ♦ a C ) - 0 (A-12) 

pp  sopvpovo 

Equation  (A-12)  evaluated  for  m^  - 0 (no  fill  gss),  becomes 

*v + <T,  - y c„  • ° <*-»> 
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Assuming  constant  Aep,  using  Eq.  (A-9),  and  noting  that,  by  definition 
of  6 


m » 6V, 
P b 
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+ C, 
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(A-14) 


(A-15) 


To  derive  an  equation  for  ur/ya,  it  is  required  only  to  manipulate 
the  expression  for  mass  average  molecular  weight  (ya  in  this  instance). 
During  burning,  propellant  gas  mixes  with  fill  gas  so  that 


a + f«V. 
o b 


u 


a 


(A-16) 


By  multiplying  Eq.  (A-16)  by  ur  and  substituting  Eq.  (A-9)  for  o^,  we 
obtain  ^ 
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The- ratios  pa/pr  and  aa/ar  are  now  given 
presented 


(A- 18) 

in  terms  of  the  ratios  already 


"JK"  W(VW“.>  <A-1« 

•J‘t  ‘ (T./T,)1/2<“r W17*  tt-20> 

In  this  report,  the  change  in  y during  burning  is  neglected,  an 
approximation  which  errs  on  the  safe  side  in  estimating  the  boob  capa- 
bility since  y actually  decreases  during  burning.  Assuming  a constant 
Ya  (equal  to  Yr ) , therefore,  results  in  a reflected  impulse  which  is  on 
the  high  side. 


In  summary,  the  equations  which  are  actually  used  in  the  calcula- 
tions are 


Pd  - (p,/145)(pa/pr) 

td  - Rt#/Z<yTf)l/2(ur/u4)1/2 


(A-21) 

(A-22) 
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*d  = R(pa/pr)(Ta/Tr)i/2(pr/pa)1/2ls/Z  (A"23) 

Z “ ci<Pa/pr)1/3/6l/3(i  “ f)1/3  (A-24) 

The  quantity,  145,  appearing  in  Eq.  (A-21)  converts  from  psi  to  MPa.  In 
the  computer  program,  values  of  ps,  ts,  and  Is  are  obtained  from  Figure  2 
by  an  interpolation  scheme  based  on  fitting  the  curves  to  a series  of 
straight  line  segments  (18  segments  for  each  correlating  curve) . 
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